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Abstract 

An analysis of the electric dipole moment (EDM) of the neutron and of the 
leptons in the minimal supersymmetric standard model (MSSM) with the most 
general allowed set of CP violating phases without generational mixing is given. 
The analysis includes the contributions from the gluino, the chargino and the neu- 
tralino exchanges to the electric dipole operator, the chromoelectic dipole operator, 
and the CP violating purely gluonic dimension six operator. It is found that the 
EDMs depend only on certain combination of the CP phases. The independent 
set of such phases is classified. The analysis of the EDMs given here provides the 
framework for the exploration of the effects of large CP violating phases on low 
energy phenomena such as the search for supersymmetry at colliders, and in the 
analyses of dark matter consistent with the experimental limits on EDMs via the 
mechanism of internal cancellations. 



It is well known that supersymmetric theories contain many new sources of 
CP violation and can produce large contributions to the electric dipole moments 
of the neutron and of the electron [1-5]. With normal size CP violating phases, 
i.e., phases 0(1), and with SUSY spectrum in the TeV range, the neutron and the 
electron EDMs already lie in excess of the current experimental limit, which for the 
neutron is|]] d n < 1.1 x 10 -25 ecm and for the electron is [[?[ d e < 4.3 x 10~ 27 ecm. 
Two approaches have usually been adopted to rectify this situation. The first is to 
make the phases small, i.e. O(10 _2 )[]I|, and the other is to use mass suppression 
by making the SUSY spectrum heavy, i.e., in the several TeV range [Q. The first 
case, however, represents fine tuning, while the second violates naturalness and 
also makes the SUSY spectrum so heavy that it may not be accesssible even 
at the LHC. Recently, a third possibility was proposed H, i.e., that of internal 
cancellations in EDMs reducing them below the experimental limits even for CP 
violating phases 0(1). 

In recent works the importance of CP violating phases on low energy phe- 
nomenona has been recognized [|| [10], 11]. In ref.|l(J it is shown that large CP 
violating phases can affect sparticle searches at colliders and in ref. [11] it is found 
that large CP violating phases can produce large effects on the neutralino relic den- 
sity consistent with the experimental constraints on the the neutron and on the 
electron EDM via the cancellation mechanism ||. However, as one goes beyond the 
framework of minimal supergravity to include non-universalities in the soft SUSY 
breaking parameters one finds that new CP violating phases arise which also affect 
low energy phenomena. Currently the effect of CP violating phases beyond the two 
CP phases allowed by the minimal supergravity cannot be investigated because 
the analytic computations of the EDMs in terms of these phases do not exist in the 
literature and consequently the EDM constraints arising from experiment cannot 
be implemented. The purpose of this Letter is to provide the analytic results for 
the EDMs beyond the minimal supergravity model by inclusion of all CP violating 
phases in the framework of MSSM. As is conventional we ignore generational mix- 
ings whose effects are known to be small. We analyse all one loop diagrams with 
the gluino, the chargino, and the neutralino exhanges for the electric dipole, and 
the chromoelectric dipole operators allowing for all phases. We also analyse the 
two loop diagrams which contribute to the purely gluonic dimension six operator 
allowing again for all CP violating phases. 

In MSSM the CP violating phases relevant for the analysis of the EDM's arise 
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from the soft SUSY breaking sector of the theory. We display this sector belowjjj]: 



V SB = mllHil 2 + m 2 2 \H 2 \ 2 - \B^e l3 E\H{ + H.c] 

+M 2 Q [u* L u L + d* L d L ] + M$u* R u R + M 2 D d R d R 
+Ml[i>:is e + ele L } + Mle* R e R 

^2m w n cos(3 1L R cos/? WL R sin/? 2yi fi J 



where (/^, q L ) are the SU(2) (slepton, squark) doublets, tan/?=| < H 2 > / < i?i > | 
where H 2 gives mass to the up quark and Hi gives mass to the down quark and 
the lepton, AVsb is the one loop contribution to the effective potential, and we 
have suppressed the generation indices. In the above, A u , A d , A e , n and B are all 
complex. Additionally, after spontaneous breaking of the electro-weak symmetry 
the vacuum expectation values of the Higgs fields are in general complex. Some 
of the phases in Eq.(l) can be eliminated by field redefinitions. However, the 
choice of which ones to eliminate is arbitrary. Rather, in our analysis we carry 
all the phases to the end and our final expressions contain only certain specific 
combinations. 

One defines the EDM of a spin-| particle by the effective lagrangian 

£i = -\d^o^F» v (2) 

In the following we will compute the contributions of the gluino, the chargino and 
the neutralino exchanges in MSSM keeping all CP violating phases. 

The gluino sector contains a phase £3 in the gluino mass term. We make a 
transformation on the gluino field to move this phase from the mass term to the 



quark-squark-gluino vertex which is then given by |L2 

: 1 — ~ ~ fc , je 3 /2^ 1 + l5~ ~ k 

i=u,d 



- £ q . 9 - 5 = V2g a T? k Y.i-^'^l^M^R + e^qlA^M^+H.c., (3) 



Here j, k — 1 — 3 are the quark and the squark color indices, a = 1 — 8 are the 
gluino color indices, and T a are the SU(3)<7 generators. The scalar fields q L and 
q R are in general linear combinations of the mass eigenstates q~i (i=l,2) so that 

q~L = Dqiiqi + D ql2 q 2 , Qr = A,21<?1 + D q22 q 2 (4) 
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where D q ij are the matrices that diagonalize the squark matrix such that D q M~ 
D q =dmg{Ml,M? 2 ), where 

M 2 _ ( M l + m i 2 + M Kl ~ Qq sin2 Ow) cos 2(3 m q (A*m - ^R q ) \ 

9 ~ ~ V m q (A q m - n*R q ) M} + m q 2 + M 2 z Q q sin 2 9 W cos 2(3 ) 

(5) ' 

Here Q u = 2/3(— 1/3) forq=u(d), R q = Vi/v^^/v*) for q=u(d), and one parametrizes 
D q so that 

/ cos | -sin|e-^\ 
^"Vsin|e^ cos % J' {b) 



where M| 21 = \M^ 2l \e %( ^ q and we choose the range of 9 q so that -y- < # g < f 

2\M? | 

where tan# g = M - 2 . In terms of the mass eigenstates q\ and g2 the gluino 

gll q22 

contribution to the EDM of the quark is given by 

<^Je - ^<)[M) - jjU^)]. (7) 



'ql J "<jl J "g2 J "<j2 
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where r* fc = e~^ 3 D q2 kD* lk , a s —^, m~ g is the gluino mass, and B(r) = (2(r — 
1) 2 ) _1 (1 + r + 2rinr{\ — r) _1 ). An explicit analysis gives T q 2 = — T^ 1 where 



M^f) = _ 9 M A m °\ A <i\ sin K - &) + |A*I sin(^ + X1 + X2 + &)\Rq\), (8) 

which holds for both signs of M| x — M| 2 , and the phases Xi (i=l,2) are defined so 
that Vi =< Hi >= \v j|e iXi (i=l,2). From Eq.(8) we see that the combinations of 
phases that enter are (a q -^) and £,3 + 9^ + xi + X2, or alternately one can choose 
them to be a q + 6» M + xi + X2 and & + M + Xi + X2- 

To discuss the contribution of the chargino exchanges we begin by exhibiting 
the chargino mass matrix 

/ \m 2 \e^ 2 V2m w sin (3e~ iX2 \ , , 

Mc ~ \ V2m w cos (3e-^ 



/i|e 



It is useful to define the transformation Mc = BrM' c B\ so that 



i\/r> _ ( l m 2| V^%smp \ , s 



v^^w sin /3 
v^m^cos/? | /i | e i ^ + « 2+ »+») 

where B R = diag(e^ 2 , e~ iXl ) and B L = diag(l, e i{ - X2+ ^). The matrix M' c can 
be diagonalized by the biunitary tranformation L r j ? M^L r L=diag (fh x +,m x +). It is 
clear that the matrix elements of Ul and Ur are functions only of the combination 
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9 — 9^ + £2 + Xi + X2- We also have U*McV~ x = diag(m x +,m x +) where U = 
(BrUr) t , and V—(BlUl)- Using the fermion-sfermion-chargino interaction we 
find that the chargino contribution to the EDM for the up quark is as follows 



— ftc., 22 m + m 2 + m 2 

€- chargmo /e = E E M^^fWiJH + (q u - g d ~)A( 1 
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(11) 



where A(r) = (2(1 — r) 2 (3 — r + 2/nr(l — r) x ) and 

= KuV^D^iU*^ - K d U* 2 D* d2k ) (12) 

and 

^"u — / — 5 e — / — v J 

\ / 2m w sin /3 v2mvi/cos/9 
Substitution of the form of U and V matrices gives: 



uil(2) 



|« u |(c OS 2 ^/2)[C/L 2i C/mJ - {+)\K^\{sin9 d )[U L2l U} i2l ]e 1 ^-^ (14) 



The terms between the brackets [ ] in Eq.(14) are functions of 9 and from the 
definition of 9 d (as given in the text following Eq.(6)) the terms between the 
brackets () in Eq.(14) are functions of the combination a d + 9^ + xi + Xi- By 
taking the imaginary part of T and using the definition of <p d (as given in the 
text following Eq.(6)) one can show that (£2 — 4>d) depends on the combinations 
(6 - a*), (6 + 0p + Xi + X2) and (a d + 9 fl + Xi + Xi)- So we are left onl y with 
the two combinations a d + 9^ + xi + X2 and £2 + + Xi + X2 with £ 2 - a d being 
just a linear combination of the first two. Similar analyses hold for the chargino 
contributions to the down quark and one gets only two phase combinations which 
are identical to the case above with a d replaced by a u . For the case of the charged 
lepton we find 

df- cha r 9 ino/e = ^ 2 m\ e £ m^M^) A(-^-) (15) 
'i/l Dill U\Y lllj ve 

where T ei = {^ e U* 2 V^) = \K e \U^ 2i ULu- A direct inspection of T ei shows that it 
depends on only one combination, i.e., ^2+^+Xi+X2- 

In order to discuss the neutralino exchange contributions we first exhibit the 
neutralino mass matrix M x o with the most general allowed set of CP violating 
phases 
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\mi\e 






|m 2 |e 



-M z sin 9 W cos /3e~ tXl M z cos W cos /3e~ lX2 
V M 2 sin 9 W sin /5e" ixi -M 2 cos 9 W sin /5e" iX2 



-M 2 sin 9 W cos /9e- ixi 
M 2 cos 9 W cos /3e" ixi 


-|/i|e ie,fl 



Next we make the transformation M x o=P^ M^ P x o where 



£2 



diag(e l 2 , e* 2 , e 



-<(^+xi) 



) 



M 2 sin ^ sin (3e~ iX2 \ 
-M z cos9 w sm(3e~ iX2 
-\n\e ie * 



(16) 



(17) 



/ 



After the transformation the matrix M^ takes the form 



\mi\ 


— M z sin9w cos/? 
V M z sin 9w sin /3e 





\rn 2 \ 
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M 2 cos cos /9e J "2" 
— M 2 cos 6V sin f3 



—M z sin 6*w cos (3 
M z cos ^vk cos (3e 1 ^ 


-Me*' 



M 2 sin 9w sin /3e 
— M 2 cos^iy sin/3 



(18) 

where 9' = + 6» M + Xi + X2, and A£ = (£1 - £ 2 )- Now the matrix M^ 

can be diagonalized by the transformation F T M^ ^=diag(m x o, m x o, m x o, m x o). 
It is clear that the transformation matrix Y is a function only of 9' and A£/2. 
Combining our results we find that the complex non hermitian and symmetric 
matrix M x o can be diagonalized using a unitary matrix X = P^oY such that 
X T M x oX=diag(m x o, m x o, m x o, m x o). We can now write down the neutralino ex- 
change contribution to the fermion EDM as follows: 



2 \ 



d E Ip 

f —neutralino I ° 



2 4 

1 ■ 2~~n 2^ 2^ Im (Vfik) 

47rsm Ow^^ 



rh y o m 2 
^fB(-AL) 



(19) 



where 



= (aoAiiDj lfc + b X 2 iD* flk + KfX bi D* f2k )(c Xi i Df 2 k - K f X bi D flk ) (20) 

where b=3(4) for T 3g = -§(§), a = -v / 2tan^(Q/ - T 3/ ), 6 = -y/2T 3f , and 
c = \/2ta,n9wQ /• We discuss now the phases that appear in the various terms 
in 7]f ik . The term proportional to a c contains the factor Xf i D*j- lk Df 2 k- It is 
easily seen that this term equals +(— )a c ^Y li (9, A£/2) sm^e"'^ 1- ^', where the 
+ (-) sign is for fc = 1(2). By doing the same analysis as for the case of the 
chargino contribution we find that the combinations that arise here are 9' , A£/2 
and aif — £i from which we can construct the three combinations: £1 + 9^ + Xi + X2, 
£2 + 9/j, +xi + X2 and a/ + 9^ +xi +X2- A similar analysis for the remaining terms of 
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Eq(20) gives exactly the same result. The sum of the gluino, the chargino and the 
neutralino exchanges discussed above gives the total contribution from the electric 
dipole operator to the quark EDM. 

The chromoelectric dipole moment d c of the quarks is defined via the effective 
dimension five operator: 

Cl = - l -d c q a ^T a qG^ a . (21) 

Contributions to d c of the quarks from the gluino, the chargino and from the 
neutralino exchange are given by 

e*»„»=^I>^)|^C(|i), (22) 

fc=l 1V1 Qk 1V1 qk 

-g 2 g s 22 fn + fn 2 + 

dq-chargino = -TfiJ- £ £ M 1 **) "^B(— §-) , (23) 
i07r fc=li=l JW <?fc iW gfe 



and 



7C _^V^ W ^ 



^-neutralino = 7^~5 J2 J2 ^(Vqik) TjT B (j7r) > ( 24 ) 
i07r fc=l i=l iW gfc iW gfe 

where B(r) is defined following eq.(7) and C(r) is given by 

^/ x 1 2r/nr 18/nr. . . 

6(r — 1) 1 — r 1 — r 

We note that all of the CP violating phases are contained in the factors Im(ri fc ), 
Im(r ? jfc), and in lm(j] q ik). But these are precisely the same factors that appear 
in the gluino, the chargino and the neutralino contributions to the electric dipole 
operator. 

Finally we look at the CP phases that enter in the CP violating purely gluonic 
dimension six operator. The gluonic dipole moment d G is defined via the effective 
dimension six operator 

C-i = —-^d G f a p 1 G aiip G P p v G 1 \ a e tJ ' 1/X ' J (26) 

where f a ^ are the Gell-Mann coefficients, ^ vXa is the totally antisymmetric tensor 
with e 0123 = +1, and G apiV is the gluon field strength. Carrying out the analysis 
including all phases we get 



d° = -Za s {^^)\m t {z\-z\)Irn{T] 2 )H{z\, 4, z^+m^-^Im^H^, z b 2 , z b )) 

(27) 
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where 

K k = e^D q2kD ; ik , 4 = (^) 2 , z q = {^f (28) 

9 9 

It is easily seen that the combination of phases involving £ 3 are similar to as for the 
gluino exchange terms discussed earlier. Similar expressions can arise if the quarks 
in the loop were from the other two generations and one gets the combinations 

The contribution to the neutron EDM using the non-relativistic SU(6) formula 
is given by d n = \(Ad d — d u ).T\ie above analysis holds at the electro- weak scale. To 
obtain the value at the hadronic scale one uses renormalization group to evolve it 
down to that scale. Thus d E = r} E d n where d E is the value at the hadronic scale and 
rj E is the QCD correction factor. The contributions of the chromoelectric dipole 
operator, and of the purely gluonic dimension six operator to the quark EDMs 
are obtained by use of the naive dimensional analysis, so that d% = -^d^r] c , and 
d G= i^ d G V G i where d% and d% are the contributions at the electro-weak scale, r/ c 
and 7]° are the QCD correction factors and M — 1.19 GeV is the chiral symmetry 
breaking scale. 

The main results of the paper are given by Eqs.(7),(ll),(15) and (19) for the 
contribution to the EDMs by the electric dipole operator, by Eqs.(22)-(24) for the 
contribution to the EDMs by the chromo-electric dipole operator and by Eq.(27) 
for the contribution to the EDM by the purely gluonic dim 6 operator. These 
formulae give the contributions to the EDMs with the most general set of CP 
violating phases with no generational mixings. 



Table 1: CP violating phases in d q and d\ in MSSM 



exchange 


u quark 


d quark 


charged leptons 


9 


OL u + 01 


otd + 0i 






6 + 


6 + 0i 




X + 


a d + 0i 


OLu + 01 






6 + 0i 


6 + 0i 


6 + 0i 


x u 


OL u + 01 


OL& + 01 


Oil + 01 




Zl + 01 


6 + 0i 


6 + 0i 




6 + 0i 


6 + 0i 


6 + 0i 


dim 6 


Oik + 01 


OL k + 0i 






k=u,d,c,s,t,b 


k=u,d,c,s,t,b 





The phases that enter in the quark and in the lepton EDM's are summarized in 
Table 1. As seen from Table 1 the electric dipole and the chromo-electric dipole 
contributions to the neutron EDMs depend on 5 phases which can be chosen to 



7 



be £j + 6*i (i=l,2,3), and + 9\ (k=u,d) where 9i=9^ +Xi+X2- The purely glu- 
onic dimension six operator contribution to the quarks depends on four additional 
phases: a:*; + Q\ (k=t,b,c,s), and thus the neutron EDM depends on nine indepen- 
dent phases. The electron EDM depends on just three independent phases: + 9\ 
(i=l,2), and a e + 9\. Thus the neutron and the electron EDM together depend 
on 10 independent phases. If we include the muon and the tau EDMs then the 
neutron and the lepton EDMs altogether depend on twelve phases, i.e., & + Q\ 
(i=l,2,3), ak + 0i (k=u,d,t,b,c,s;e,/i, r). If we retain only the dominant top-stop 
contribution to the purely gluonic dimension six operator, then the total number 
of phases reduces from 12 to 9. 

For the case of minimal supergravity, all evolve from the phase of A Q 
where A Q is the common value of Aj at the GUT scale. Similarly all & evolve from 
the phase £i of the universal gaugino mass mi at the GUT scale. In this case the 

^2 2 

ten phase combinations that appear in d n and d e collapse to just two independent 

ones: ao+9\, and £i+0i. Often one sets £i = 0, Xi — = X2 and chooses the 

2 2 

independent phases in minimal supergravity to be cto and 9^. In that limit the 
results of this Letter limit to the results of ref.||. 

The analysis presented here is the first complete analysis of the neutron and 
of the lepton EDMs with all allowed CP phases in MSSM under the restriction of 
no generational mixings. In ref. consistency with large CP violating phases was 
achieved by internal cancellations to satisfy the experimental EMD constraints 
with just two CP phases. In the present analysis since the electron EDM depends 
on three independent phases and the neutron EDM depends on nine independent 
phases, the satisfaction of the experimental EDM constraints with large CP vio- 
lating phases can be achieved over a much large region. The analysis provides the 
framework for the investigation of the effects of large CP violating phases on low 
energy physics. 

This research was supported in part by NSF grant PHY-96020274. 
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